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ABSTRACT 


Problems  of  transmission  and  reflection  of  long,  time  harmonic,  free 
surface  gravity  waves  obliquely  incident  from  a  constant  depth  ocean  upon 
linearly  varying  bottom  topographies  are  considered.  A  solution  to  the 
vertically  integrated  dynamic  equations  over  linear  depth  variation  is 
developed  in  terms  of  Kummer  functions.  Coriolis  effects  are  included  but 
primary  interest  is  on  Class  I  (high  frequency)  long  waves.  Three  cases  are 
treated;  the  continental  slope  and  shelf,  the  submerged  ridge  and  the 
submerged  trench.  \ 
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Introduction 


This  study  considers  the  effects  of  three  bottom  topographies;  the 
continental  shelf  and  slope,  submerged  ridges  and  submerged  trenches,  on 
originally  plane  propagating  long  free  surface  gravity  waves.  In  each  case, 
linearized  long  wave  theory  for  a  homogeneous  perfect  fluid  is  used  on  one 
dimensional  bottom  topographies  (parallel  contours)  which  possess  constant 
and  linear  cross  sectional  depth  variations.  The  Coriolis  effect  due  to  the 
earth's  rotation  has  been  included  in  the  theory.  A  distinction  may  be  made 
on  the  basis  of  the  wave  frequency,  u>,  and  the  Coriolis  frequency,  f,  between 
the  two  classes  of  long  water  waves.  The  Class  I  or  inertiogravitational 
waves  are  characterized  by  m  »  f.  These  waves  can  propagate  in  both 
directions  along  a  N-S  coast  albeit  with  differing  phase  speeds  for  non-zero 
rotation  and  have  periods  which  fall  in  the  tsunami  range  and  shorter,  i.e., 
on  the  order  of  an  hour  or  less.  The  Class  II  or  quasi -geostrophic  waves  are 
characterized  by  u  «  f.  These  waves  propagate  in  only  one  direction  and 
cannot  exist  in  the  absence  of  rotation.  They  are  an  oceanographic  scale 
phenomena  with  periods  on  the  order  of  many  hours  or  even  days.  This  study 
is  directed  at  the  Class  I  waves  where  the  Coriolis  effect  plays  only  a 
modifying  role.  Leblond  and  Mysak  (1978)  provide  an  important  overview  of 
the  general  area  of  ocean  waves. 

The  solution  utilized  over  linear  depth  variation  was  noted  by  Shaw 
(1974)  and  Guza  and  Davis  (1974)  in  terms  of  Kummer  functions,  but  no 
numerical  results  were  given  until  Shaw  (1977)  due  to  an  apparent  lack  of 
tables  and/or  subroutines  for  these  functions.  The  application  of  this 
linear  topography  solution  considered  is  that  of  long  wave  scattering  by 
topographical  features.  Since  these  problems  assume  a  progressive  plane 
wave  field  incident  from  a  semi-infinite  ocean  onto  the  feature  to  be  studied. 
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the  wave  may  not  be  trapped  (within  a  linear  theory).  Wave  trapping  is  studied 
in  a  related  paper,  Shaw  and  Neu  (1981).  Of  interest  here  is  the  proportion 
of  the  incident  wave  energy  reflected  by  or  transmitted  through  the  feature  and 
the  resulting  wave  amplitudes,  especially  at  a  coastline  for  the  shelf-slope 
case. 

The  case  of  an  infinite  shelf,  i.e.,  a  slope  connecting  two  constant  depth 
regions  with  no  coastline,  was  carried  out  by  asymptotic  methods  by  Voyt  and 
Sebekin  (1974)  for  an  arbitrary  monotonic  depth  in  the  slope  region  with  the 
particular  case  of  normal  incidence  and  a  linear  slope  carried  out  analytically. 
This  problem  has  also  been  studied  by  Dean  (1964)  where  the  analytical 
solution  for  a  linear  slope  and  normal  incidence  was  given  without  Coriolis 
effects  and  coastline.  Shaw  (1977)  considered  the  case  of  long  waves  obliquely 
incident  upon  a  linear  slope  and  constant  depth  shelf  which  is  terminated  by  a 
partially  reflecting  vertical  wall.  Recently,  Abe  and  Ishii  (1980)  presented 
a  similar  solution  to  oblique  incidence  on  a  linear  slope  and,  using  a  constant 
depth  shelf,  treat  transmission  and  reflection  problems  in  the  cases  of  an 
infinite  shelf  and  a  shelf  terminated  by  a  rigid  vertical  wall.  Refraction, 
diffraction  and  scattering  are  treated  from  a  ray  theory  viewpoint  by  Murty 
(1977),  who  also  provides  a  review  of  related  work. 

This  study  will  develop  the  basic  solution  for  long  waves  over  a  linearly 
varying  bottom  topography  and  apply  this  solution  to  problems  of  the  transmission 
and  reflection  of  waves  obliquely  incident  upon  submerged  ridges  and  trenches 
as  well  as  the  continental  slope  and  shelf  where  both  slope  and  shelf  have  linear 
slopes  and  the  shelf  terminates  with  zero  depth  at  the  coastline. 
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Basic  Formulation 

The  topographical  features  to  be  studied  are  assumed  to  have  contours 
parallet  to  a  y  axis  and  a  piecewise  continuous  cross  section  composed  of 
segments  whose  depth  varies  linearly  in  an  x  direction.  Although  rotation  on 
an  f  plane  is  included,  the  emphasis  will  be  on  first  class  surface  gravity 
waves  where  Coriolis  effects  are  a  modifying  rather  than  fundamental  influence. 

The  vertically  integrated  equations  governing  long  wave  motion  in  water 
of  variable  depth,  H,  and  Coriolis  parameter  f  are,  Lonquet-Higgins  (1968), 


0 


(4) 
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which  will  have  dii  erent  solutions  for  different  topographies,  F(>T). 

In  each  of  the  cases  to  be  considered,  it  is  assumed  that  the  topographical 
feature  is  preceded  by  a  region  of  constant  depth  over  which  the  incident 
wave  field  travels.  The  depth  H-j  shall  be  used  as  the  reference  depth.  Thus 
in  this  region,  h(x)  =  1  and  equation  (4)  becomes, 

V2  ?  +  o'  ?  =  0  (5) 

o  V2 

where  we  have  set  a  =  u  [1  -  7C]  .  A  plane,  unit,  incident  wave  field  whose 

normal  makes  an  angle  e  with  the  x  axis  is  assumed  of  the  form. 

If.  =  exp  [i  a  cos 0  x  -  i  a  sino  yl. 

Also  present  in  this  region  will  be  a  corresponding  wave  reflected  from  the 
topographical  feature,  traveling  in  the  negative  x  direction  and  given  by 

?r  *  R  exp[- i  a  cose  x  -  i  a  sino  y] 

where  R  is  the  coefficient  of  reflection.  Thus,  the  total  wave  field  in  this 
region  is, 

jf  =  [exp(i  a  cose  x]  +  R  exp(-i  a  cose  x)  ]  exp(-i  a  sine  y)  (6) 

which  satisfies  equation  (5).  It  is  noted  here  that  solutions  in  any  region 
must  have  the  same  y  dependence,  therefore  it  will  be  dropped  from  future 
equations  along  with  the  explicit  harmonic  time  dependence  which  has  already 
been  dropped. 

In  general,  for  a  region  of  constant  non-dimensional  depth,  hc,  the 
governing  equation  beccres, 

—  2  —  ?  — 
hc  7C  c  +  C  = 


0 


(7) 
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i  2  1/2  _ 

+  C0  exp[i  a( —  -  sin  e)  x] 
c  h 

c 

where  and  C£  are  arbitrary  constants. 

In  regions  of  linearly  varying  non-dimensional  depth  defined  by 

h(x)  =  5  +  yx  (9) 


where  6  and  y  are  constants  {y  f  0),  equation  (4)  becomes. 


(6  +  yTj  V2  Z  +  y  f§  +  i  Ty  +  a2  5  =  0 


Introducing  a  new  variable,  z  =  5  +  yx  and  substituting  for  the  y  derivatives 
of  Ti  this  equation  is , 


2—  —  2  —  2  2 
d  5  ^  dq  ,  ra  +  f  y  a  sine  -  a  sin  e  z,  - 

~TZ  Hz  +  [ - - - 1  t 

dz  y 


01) 


which  is  reducible  to  the  confluent  hypergeometric  equation  as  seen  in 
Chapter  6,  section  2  of  Erdelyi,  et  al .  (1953).  Two  distinct  cases  must  be 
considered.  If  e  =  0,  equation  (11)  corresponds  to  equation  (8)  of  this 
reference  and  the  solutions  are, 

CQ  (2aZ1/2/yl 

where  CQ  represents  any  two  independent  solutions  to  the  zero  order  Bessel 
equation.  If  5  f  G,  v;e  have  the  form  of  equation  (6)  of  this  reference  with 
solutions , 


exp  [-  ?/2]  .  [F  a  ,c,c  ] 
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where, 

_  y  sine  -  a  -  T  y  sine  _  1  -  7  a 

2y  sine  2  2y  sine 

c  =  1 

_  2a  sine 
5  y 

F  represents  any  two  independent  solutions  to  the  confluent  hypergeometric 
equation.  Since  c  =  1,  one  of  the  solutions  will  be  of  logarithmic  form. 

This  form  is  appropriate  for  c  >  0.  For  5  <  0,  the  Kummer  transformation  leads 
to  the  form, 

exp[+  c/2]  .  F[a' ,  1,-5] 

where  s’  ■  1  -  a  -  w  • 

The  solution  in  regions  of  non-zero,  constant  slope  bottom  topography, 
for  e  /  0,  is  therefore  (using  the  notation  of  Abramowi tz  and  Stegun,  1964), 

1  =  expt-c/2]  {c3M(a,l,c)  +  c4U(a,l,c)};  c  >  0 

(12) 

=  exp[+  c/2]  tc^Ka1 ,1  ,-c)  +  c4U(a,l  ,-c) } ;  C<  0 

where  c^  and  c4  are  arbitrary  constants.  M  and  U  represent  the  two  independent 
Kummer  functions.  Programs  to  evaluate  these  functions  have  been  written  by 
and  are  available  from  the  authors.  These  solutions  may  also  be  expressed  in 
terms  of  Whittaker  functions. 


r  ' — - — 
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1 

Continental  Slo 

pe  and  Shelf  I 

The  first  application  considered  is  the  reflection  of  a  unit  plane  wave 
from  a  continental  slope  extending  from  x  =  0  to  x  =  A  and  sloping  shelf 
extending  from  x  =  A  to  B.  The  shelf  terminates  with  zero  depth  at  the  coastline, 
x  =  B.  A  cross  section  of  the  ocean,  shown  in  figure  1,  is  divided  into  three 
regions , 

Region  I:  -»  <  x  <  0;  H(x)  =  H-j 

Region  II:  0  <  x  <  A;  H(x)  =  H-]  -  [H-j  -  H2)  x/A 

Region  III:  A  <  x  <  B;  H(x)  = 

The  reference  length  is  chosen  to  be  A  and  h 2  =  H2/H-j ,  b  =  B/A  are  defined. 

The  solution  in  region  I,  as  given  by  equation  (6)  is, 

Cj  =  exp(i  a  cose  x")  +  R  exp(-i  a  cose  x) ;  <  x  <  0  (13) 

In  region  II,  ¥(7)  =  z  =  1  -  (1  -  F2)  x.  Therefore,  in  equation  (9), 
y  =  -  (1  -  F2)  and  the  solution  in  this  region  (e  t  0),  from  equation  (12),  is, 

Fjj  =  exp[-£2/2]  {A2H(a2,l ,c2)  +  B2U(a2,l ,c2) } ;  e  <  0 

0  <  x  <  1  (14) 

=  exp[s2/2]  {A2M(a£,l,-s2)  +  B2U(a£  ,1  ,-c2) } ;  e  >  0 
where  A2  and  B2  are  constants  and, 


a  -  1  -  T  , 
a2  2 


a>  =  1  _  a  -  L±_L 
a2  a2  2 


2(1  -  h2)sine 

1  +  f _ a 

2(1  -  h2)sine 


=  2a  sine  [x - — ]  . 

£  1  =  h2 
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In  region  III,  h(x)  =  z  =  h9  ,  y  = - - — ,  and 

2  F  -  1  F  -  1 


’III  =  exP[--3/2J  {A3M(a3,1 ,s3)  +  B3U(a3,l  ,c3)} 


=  exp[?3/2]  {A3M{a^,IJ-c3)  +  B3U(a^,l,-c3)} 


e  <  0 

1  <  7  <  F 
6  >  0 


(15) 


where  A3  and  E>3  are  constants  and, 
-  1  -  I  , .  a(F-  1) 


3  2 


2h3  sine 


a '  =  l  - 

3  3  2  2F2  sine 


C3  =  2a  sine  [x  -  b] 


However,  as  x  -*•  F  (^3  -*■  0) ,  U  ->■  -»  .  Therefore,  to  insure  a  bounded  solution  at 
the  coastline,  B3  must  be  set  to  zero.  The  solution  in  region  III  is  then, 

tni  =  exp[-  53/2]  A3M(a3,l,c3);  e  <  0,  1  <  7  <  (16) 


=  exp[c3/2]  A3M(a£,  1  -  £3) ;  0  >  0 

The  evaluation  of  the  coefficients  R,  Ag,  and  A,  is  accomplished  by 
requiring  continuity  of  the  free  surface  elevation  and  mass  flux  at  the  inter¬ 
faces  x  =  0  and  x  =  1.  This  introduces  the  following  matching  conditions, 

(x  =  0)  =  Fn(x  =  0) 

3^1  _  3Ft7  _ 

—  (x  =  0)  =  U  =  0) 

3X  ;  X 

tjj(x  1)  tjjj(x  -  1) 


V 


(17) 
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3? 


II  #- 


(x  =  D  = 


3? 


Ill 


(X  =  1) 


3X 


3x 


Equations  (17)  are  four  simultaneous  algebraic  equations  in  R,  A^,  B^,  and 
A^  which  may  be  solved  to  complete  the  solution  to  the  continental  slope  and 
shelf  problem. 


V 
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The  second  problem  considered  is  that  of  the  transmission  and  reflection 
of  a  unit  plane  wave  by  a  submerged  ridge,  which  we  divide  into  four  regions: 

Region  I:  -»  <  x  <  -  A;  H(x)  =  H-j 

Region  II:  -A  <  x  <  0;  H(x)  =  H2  -  [H^  -  x/A 

Region  III:  0  <  x  <  B;  H(x)  =  H2  +  [H^  -  Hg ]  x/B 

Region  IV:  B  <  x  <  +®;  H(x)  =  H3 

as  shown  in  Figure  2.  Again,  A  is  used  to  scale  horizontal  distance  and, 

F2  =  H2/Hi  ,  F3  =  H3/H-J ,  and  F  =  B/A  are  defined. 

In  region  I,  the  solution  is  again  given  by  equation  (13)  but  applies  to 

-a.  <  x  <  -1.  In  region  II,  F(x)  =  F2  -  (1  -  F2)  F,  y  =  -  (1  -  F2)  and  the 

solution  is 


=  exp[-c2/2]  {A2M(a2,l  ,c2)  +  B2'J(a2 , 1  ,f;2) >  *>  0  <  0  ,  -1  <  x  <  0  (18) 

=  expU2/2]  {A2M(a2,l  ,-C2)  +  B2U(a2,l  ,-c2)>  ;  0  >  0 

where. 


_ a _ 

2(1  -  F2)  sine 


g.  _  1  +  f _ a _ 

2  2  2(1  -  F2)  sina 

-  ^2 

-  2a  sine  (x - ] 

i  1  -  h2 

In  region  III,  F(x)  =  F.  +  (F^  -  F^)  F/F,  y  =  (F^  -  F2)/F  and  the  solution  is 


C 1 1 1  =  exp [ C3/2 ]  {A3M(a^,l ,-C3)  +  B3U( a^ , 1 »-C3) )  ;  8<0,  0<x<F  (19) 

=  exp[-£3/2]  {A3M(a3,l ,^3)  +  B3U( a3 ,1  ,c3) )  ;  0  >  0 


f 


where , 
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1  -  7 


a3  2 


2(h0  -  h0)  sino 

o  L 


i '  =  i-J-t  + - 


2(F^  -  F^)  sine 


5-  =  2a  sine  [x  +  = — ~r] 
h3  -  h2 


In  region  IV,  h(x)  =  h^,  a  constant.  If  only  the  right-running  transmitted 
wave  is  allowed  in  this  region,  the  solution  is  given  by 

1  2  1/2-  -  - 

Civ  =  A4  exPt'1"  a  ( - sin  e)  x]  ;  b  <  x  <  ■»  (20) 

^3 

1  2  1 

If  H_  >  H, ,  then  —  <  1  and  if  sin  e  >  —  ,  this  solution  becomes  a  decaying 
h3  h3 

exponential.  In  this  case,  no  energy  is  transmitted  in  region  IV.  The  wave 

is  totally  reflected  by  the  ridge.  The  critical  incident  angle,  ecr>  for  such 

total  reflection  is  defined  by  sine  =  ~~  ,  which  agrees  with  Snell's  law. 

*3 

Thus,  the  solution  in  region  IV  is, 

—  12  —  1  2  _  _ 

Civ  =  A4  exPtia  -  sin  e)  X]  —  >  sin  e  i  <  x  <  «  (21) 

h3  *3 


2  1  */ c  _  I  2 

=  Aa  exp[-a  (sin  e - )  x)  —  <  sin  e 

F3  h3 

To  determine  the  coefficients  R,  A^,  B^,  A^,  and  A^,  continuity  of  free 
surface  elevation  and  mass  flux  across  the  interfaces  is  once  again  required. 


f 
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Submerged  Trench 


The  submerged  trench  problem,  Figure  3,  is  very  similar  to  the  submerged 
ridge  problem.  The  four  regions  are  defined  by 

Region  I:  —  <  x  <  -A;  H(x)  =  Hj 

Region  II:  -A  <  x  <  0;  H(x)  =  H2  +  [H2  -  H-j]  x/A 

Region  III:  0  <  x  <  B;  H(x)  =  H2  -  [H2  -  H3]  x/B 

Region  IV:  B  <  x  <  +  H(x)  =  H3 

Again,  A  is  used  to  scale  horizontal  distance  and  K2  =  H2/H^ ,  Tf3  =  H3/H^  and 
b  =  B/A  art  defined.  The  solutions  in  each  region  are 


Cj  =  exp(i  a  cose  x)  +  R  exp(-i  a  cose  x)  ;  -<*>  <  x  <  -1 


(22) 


CH  =  exp[c2/2]  {A2M(a^,l,-c2)  +  B2U(a‘ ,1 ,-s2) }  ;  e  <  0  ,  -1  <  x  <  0  (23) 

*  exp[-c2/2]  (A2M(a2,l,s2)  +  B2U(a2,l ,c2) }  ;  e  >  0 


a  =  1  -  f 
a2  2 


— — -  ,  ai  =  1  -  a? 

2(h,  -  1)  sine  c  1 


So  =  2  a  sine  [x  + 


f2  -  1 


CIII  =  exP("?3^1  {A3M(a3,l,c3)  +  B3U(a3,l,s3)}  ;  6<0,  0<x<b  (24) 

=  exp[s3/2]  {A3M(a^,l,-e3)  +  B3U(a^,l ,-c3)l  ;  e  >  0 


a3  = 


-I— I— L.  +  _ 2^ - - -  a*  =i  _  a 

2(h2  -  h3)  sine  J  J 


F?F 

S-  =  2a  sine  [x  -  — — j 
h2  -  h3 
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Results 


Continental  Slope  and  Shelf 


The  effect  of  rotation  is  considered  first.  Figure  4  is  a  plot  of  the 
amplification  of  the  incident  wave  at  the  coast  versus  the  incident  angle,  e, 
for  periods  (T)  of  900,  1200,  2400,  and  4800  seconds,  each  over  the  same  topography 
at  latitudes  of  0  and  60  degrees.  At  zero  latitude,  where  no  Coriolis  effect  is 
present,  the  curves  are  symnetric  about  0=0.  As  rotation  is  introduced,  the 
curves  shift  slightly  and  become  asymmetric.  It  is  interesting  to  note  that 
the  2400  second  period  curve  shifts  in  a  direction  opposite  that  of  the  other 
three.  It  is  seen  that  the  Coriolis  effect  is  larger  for  the  longer  period 
waves  than  for  the  shorter;  however  even  at  the  rather  large  latitude  of  50 
degrees,  the  effect  is  small  for  the  class  I  type  wave  periods  of  interest  here. 

All  further  results  therefore  will  be  taken  without  Coriolis  effects,  i.e.,  f 
equal  to  zero. 

It  is  now  convenient  to  define 

X  =  T  ✓gHJ  (26) 

where  T  is  the  wave  period.  This  X  is  the  wavelength  at  the  edge  of  the 
shelf  (x  =  A).  It  was  found  that  the  amplification  at  the  coast,  as  a  function 
of  incident  angle,  is  dependent  upon  the  ratio  of  the  shelf  length  to  x, 

B  -  A  — 

— - —  ,  and  the  ratio  hg.  There  is  also  a  slight  dependence  upon  the  length  A. 

The  amplification  at  the  coast  at  zero  incident  angle  is  plotted  versus  — ^ — 

for  values  of  l/hg  of  15,  25,  and  35  in  figure  5.  This  curve  clearly  shows  the 

B  -  A 

cyclic  behavior  and  resonance  peaks  of  the  amplifications.  As  ~ —  is  increased, 

the  cycle  repeats  at  intervals  of  just  under  0.25  with  each  cycle  slightly  higher 
and  with  a  larger  peak  than  the  last.  The  amplitude  of  the  peaks  also  increases 
with  l/F^. 
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Figures  6  and  7  illustrate  the  angular  dependence  of  the  amplification. 

They  are  plots  of  the  amplification  at  the  coast  versus  the  incident  angle 

through  one  cycle  of  Figure  5,  from  — j —  =  0.8  to  1.05  in  steps  of  0.025  and  at 

1/^2  =  25.  Only  6  _<  0  is  shown  due  to  the  symmetry.  The  general  shape  of  the 

8  -  A 

curves  repeat  with  each  cycle.  The  curves  from  — - —  =  0.8  to  0.925,  shown  in 

Figure  6,  are  monotonic  with  the  maximum  at  9  =  0.  These  curves  correspond  to 

points  to  the  left  of  a  resonance  peak  in  Figure  5,  on  an  upward  slope.  In 
Figure  7,  the  curves  correspond  to  points  to  the  right  of  a  resonance  peak  on 
a  downward  slope.  They  are  no  longer  monotonic  in  e,  but  have  maximums  away 

from  zero  near  e  *  -15,  -45,  -60,  and  -75  for  ^  =  0.95,  0.975,  1.0,  and  1.025 

respectively.  The  curve  — ^ =  1.05  corresponds  to  the  — - —  =  0.8  curve  and 
completes  the  cycle. 

The  differences  in  the  slopes  of  the  curves  to  the  left  and  to  the  right 

of  a  resonance  peak  may  be  heuristically  explained  as  follows.  As  the  magnitude 

of  0  increases,  the  wavelength  perpendicular  to  the  coastline,  corresponding  to 

A  in  Figure  5,  also  increases.  Thus,  the  ratio  of  the  shelf  length  to  the 

B  -  A 

perpendicular  wavelength  is  always  less  than  or  equal  to  — - —  .  Referring  to 

B  -  A 

Figure  5,  at  a  point  to  the  left  of  a  resonance  peak,  e.g.,  — ^ —  =  0.9,  the  ratio 
of  shelf  length  to  perpendicular  wavelength,  when  j  e  |  >0,  is  to  the  left  of  0.9  in 
a  region  of  less  amplification.  At  a  point  to  the  right  of  the  peak,  e.g., 

— — =  0.975,  the  shelf  length  to  perpendicular  wavelength  is  in  a  region  of 

A 

greater  amplification.  These  waves  encounter  a  resonance  peak  at  an  angle 

B  -  A 

away  from  zero.  In  the  case  of  — - —  =  0.975,  it  is  near  9  =  -45.  These 

A 

maximum  amplifications  away  from  6=0  are  more  predominant  in  waves  of  shorter 
periods  where  the  resonance  peaks  of  Figure  5  are  larger.  Of  course,  this  so 
called  perpendicular  wavelength  is  difficult  to  define  due  to  the  fact  that  the 
wave  is  turning  in  toward  the  coast  as  it  moves  up  the  slope  and  shelf,  however 
a  qualitative  understanding  of  the  phenomenon  may  be  reached  through  the  above 
argument. 
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Submerged  Ridge 

Figures  8  and  9  are  plots  of  the  magnitude  of  the  reflection  coefficient,  R, 
versus  the  incident  angle  for  several  periods.  The  transmission  coefficient  is 
inferred  by  R  since  energy  must  be  conserved.  The  ridge  is  symmetric  with 
dimensions  roughly  coinciding  with  those  of  the  Mid-Atlantic  Ridge;  H-j  =  H3  *  5  km, 
Hn,  -  2  km,  A  =  B  ■  100  km.  For  waves  whose  lengths  are  on  the  order  of  the  total 
width  of  the  ridge,  the  behavior  is  similar  to  the  T  =  900  sec.  curve.  As  the 
period  (and  wavelength)  is  increased,  a  point  of  zero  reflection  develops  and  the 
ridge  acts  similar  to  a  band-pass  filter.  This  point  of  zero  reflection  moves 
inward  along  the  incident  angle  axis  with  increasing  period  until  reaching  a 
point,  ]  e  | 53°  in  this  case,  where  it  becomes  relatively  stationary.  The  shape 
of  the  curves  in  the  region  where  |  e  |  is  less  than  the  point  of  zero  reflection 
changes  steadily  from  a  shape  of  which  the  T  =  1200  sec.  curve  is  typical,  where 
the  reflection  rises  to  a  local  maximum  then  drops  off  as  e  approaches  0,  to  a 
monotonically  increasing  curve  as  in  T  =  2100  sec.  As  the  period  is  increased 
beyond  this  point,  the  shape  of  the  curves  remains  fairly  constant  with  the 
reflection  at  0  =  0  rising  to  a  maximum  then  slowly  falling  off  as  seen  in  Figure  9. 

In  Figure  10,  an  asymmetric  ridge  is  considered,  illustrating  the  possibility 
of  total  reflection.  Here,  =  5  km,  =  2  km,  =  6.5  km,  A  =  100  km,  and 
B  =  120  km.  The  critical  angle  for  an  incident  wave  traveling  in  region  I  to  be 
totally  reflected  is  61.29°.  For  angles  of  magnitude  less  than  the  critical  angle, 
the  curves  retain  a  shape  similar  to  the  corresponding  curves  in  the  symmetric  case. 
Notice  that  the  rise  to  total  reflection  is  very  steep,  most  of  which  occurs  in  an 
interval  of  about  10  degrees. 


I 
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Submerged  Trench 

Again,  the  magnitude  of  the  reflection  coefficient  is  plotted  versus  the 
incident  angle  for  several  values  of  T  in  Figures  11  and  12.  The  trench  is 
symmetric  with  =  5  km,  Hg  =  9  km,  A  =  B  *  100  km.  As  in  the  ridge 

problem,  there  is  a  point  of  zero  reflection,  however  it  is  present  at  the 
period  of  900  sec.  occurring  at  about  |  e  |  =  35°.  As  the  period  is  increased, 
the  point  of  zero  reflection  moves  out  to  about  |  e  }  =  39°  where  it  remains 
relatively  stationary.  Notice  also  that  the  T  =  900  sec.  curve  has  very  little 
reflection  for  incident  angles  of  magnitude  less  than  the  point  of  zero 
reflection.  The  other  major  difference  from  the  ridge  problem  is  that  at  sharp 
angles,  the  smaller  period  waves  are  almost  totally  reflected.  This  is  due  to 
the  fact  that  they  are  turned  back  out  away  from  the  trench  by  the  slope  of 
region  II.  As  the  ratio  of  Hg  to  H-|  is  increased,  this  reflection  will  increase. 

In  Figure  13,  an  asymmetric  trench  is  considered.  Although  the  possibility 
of  total  reflection  exists  in  the  trench  problem  when  >  H-|  as  in  the  ridge 
problem,  a  model  of  the  Japan  Trench  was  chosen  here  with  <  H^.  The  dimensions 
are:  =  6  km,  =  9.5  km,  =  3  km,  A  =  60  km,  B  =  90  km.  Note  that  in  this 

case,  only  the  T  =  1800  sec.  curve  has  a  point  of  zero  reflection.  The  other 
curves  turn  back  up  before  reaching  zero.  In  general,  the  reflections  are  larger 
than  in  Figures  11  and  12  for  the  angles  of  lesser  magnitude,  due  to  the  shallower 
depth  H^,  and  smaller  for  angles  of  greater  magnitude  due  to  the  smaller  ratio 
of  H2  to  H-j . 


Conclusion 


This  presentation  deals  with  an  analytical  study  of  a  set  of  physical 
problems.  No  comparison  is  made  at  this  time  to  observed  oceanographic  behavior. 
Long  wave  measurements  across  ridges  and  trenches  are  sparse;  however  it  is  clear 
that  long  waves  actually  do  traverse  ridges  and  trenches  as  they  cross  the  ocean, 
e.g.,  tsunamis,  and  some  reflection  and  transmission  of  energy  must  occur. 

The  amplification  of  wave  heights  at  a  coastline  is  more  readily  observed, 
but  deep  water  wave  amplitudes,  measurements  from  which  these  coastal  values  are 
obtained,  are  again  sparse  --  even  for  the  important  case  of  tsunamis.  The  present 
analytical  solution  presents  a  "transfer  function"  representing  the  effect  of  the 
continental  shelf-slope  topography  on  an  incident  wave,  i.e.,  the  effect  of  shelf 
resonances  by  selectively  amplifying  certain  frequencies.  This  transfer  function 
transforms  the  deep  water  spectrum  to  a  coastline  spectrum  for  a  plane  incident 
wave  but  is  dependent  on  the  incident  angle  as  well  as  the  topography.  Assuming 
the  topography  to  be  known,  knowledge  of  the  coastline  spectrum  and  the  original 
incident  angle  would  allow  reconstruction  of  the  deep  water  spectrum.  Abe  (1981) 
has  attempted  to  identify  an  incident  angle  but  for  tsunami  sources  on  the  same 
shelf  as  the  observation  points;  the  reconstruction  here  would  require  a  distant 
source  mechanism.  Such  problems  are  presently  being  studied. 
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